Josephson effect through a multilevel dot near a singlet-triplet transition 



Minchul Lee,' Thibaut Jonckheere,^ and Thierry Martin^' 

'Department of Applied Physics, College of Applied Science, Kyung Hee University, Yongin 449-701, Korea 
^Centre de Physique Theorique, UMR6207, Case 907, Luminy, 13288 Marseille Cedex 9, France 
^Universite de la Mediterranee, 13288 Marseille Cedex 9, France 

We investigate the Josephson effect through a two-level quantum dot with an exchange coupling between two 
dot electrons. We compute the superconducting phase relationship and construct the phase diagram in the super- 
conducting gap-exchange coupling plane in the regime of the singlet-triplet transition driven by the exchange 
coupling. In our study two configurations for the dot-lead coupling are considered: one where effectively only 
one channel couples to the dot, and the other where the two dot orbitals have opposite parities. Perturbative 
analysis in the weak-coupling limit reveals that the system experiences transitions from to tt (negative critical 
current) behavior, depending on the parity of the orbitals and the spin correlation between dot electrons. The 
strong coupling regime is tackled with the numerical renormalization group method, which first characterizes 
the Kondo correlations due to the dot-lead coupling and the exchange coupling in the absence of superconduc- 
tivity. In the presence of superconductivity, many-body correlations such as two-stage Kondo effect compete 
with the superconductivity and the comparison between the gap and the relevant Kondo temperature scales al- 
lows to predict a rich variety of phase diagrams for the ground state of the system and for the Josephson current. 
Numerical calculations predicts that our system can exhibit Kondo-driven O-vr-O or tt-O-tt double transitions and, 
more interestingly, that if proper conditions are met a Kondo-assisted tt -junction can arise, which is contrary to 
a common belief that the Kondo effect opens a resonant level and makes the 0-junction. Our predictions could 
be probed experimentally for a buckminster fullerene sandwiched between two superconductors. 

PACS numbers: 73.63.-b, 74.50.+r, 72.15.Qm, 73.63.Kv 



I. INTRODUCTION 

The Josephson effect''^ is one of the most celebrated man- 
ifestation of many body correlations in condensed matter 
physics: a Cooper pair cuiTent^ between two bulk supercon- 
ductors separated by an intermediate region or arbitrary na- 
ture can flow even in the absence of an applied bias. Over 
the last few decades, the Josephson effect has become a very 
active field of theoretical"*"'^ and experimental"^^- investiga- 
tion in the context of mesoscopic devices, devices which are 
small enough that electron transport occurs in a phase coher- 
ent manner. Because the tunneling of Cooper pairs through 
the junction is greatly affected by the physical properties of 
the segment between superconducting electrodes, the study of 
the Josephson current provides a new way to investigate the 
electronic properties of the medium. In early days, thin lay- 
ers of insulators and metals were used to form the Josephson 
junctioni^ Advance in nanofabrication technology can now 
enable one to make the middle segment small enough to be 
considered as a quantum dot (QD), a zero dimensional entity 
bridging the two superconductorsi^^^ Furthermore, even 
a real (or artificial) molecule can be inserted between two 
closely positioned superconducting leads to form a molecular 
Josephson junction (MJJ)ji^i2i^ Quantum dots connected to 
normal metal leads are known to exhibit the Coulomb block- 
ade phenomenon due to their large charging energies. In- 
terestingly, at low enough temperatures, when an odd num- 
ber of electrons occupy the QD, one can reach the Kondo 
regimei^^ Yet the leads can be chosen to be superconduc- 
tors which lead to a competition between Kondo physics and 
superconductivity. '''"'^•^'•^^ The purpose of the present work 
is precisely to study the Josephson effect through a multilevel 
quantum dot in this context, with applications to molecular 



spintronics. 

Indeed, the QD-JJ have received a great theoretical and ex- 
perimental attention because they can exhibit an interesting 
competition between two many-body correlations: the super- 
conductivity and the Kondo effect. Due to its small size, the 
QD has a large Coulomb charging energy. The Kondo ef- 
fect then emerges for such a small QD coupled strongly to the 
leads when the QD has a localized magnetic moment, that is, 
nonzero total spin of electrons in it. At temperatures below the 
so-called Kondo temperature Tk?^ the conduction electrons 
in the leads screen the localized moment through multiple co- 
tunneling spin-flip processes, forming a spin-singlet ground 
state, and induces a resonance level at the Fermi energy, which 
increases the linear conductance up to the unitary-limit value 
{— jh) that is otherwise completely suppressed due to 
the strong Coulomb repulsion. If the leads consist of s-wave 
superconductors, the conduction electrons form spin-singlet 
Cooper pairs incapable of flipping the QD spin. It has been 
known,"*"^'^" and recently probed,^''^^'^^ that in the weak dot- 
lead coupling limit the large Coulomb repulsion only allows 
the electrons in a Cooper pair to tunnel one by one via vir- 
tual processes in which the spin ordering of the pair is re- 
versed, leading to a tt junction, and that the localized moment 
remains unscreened. In the opposite limit where the Kondo 
temperature exceeds the superconducting gap A, however, the 
induced Kondo resonance level restores the junction state of 
the supercuiTent.^''"'""' As a result, one can drive a phase tran- 
sition between spin singlet (0 junction) and doublet (tt junc- 
tion) states by changing the relative strengths of Tk and A. 

Current issues about electronic transport through a QD or 
a molecule go beyond the spin-degenerate single-level model 
and take into account multi-level structures and/or possible 
magnetic interactions. For example, the theoretical predic- 
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tion that the two-level quantum dots (TLQDs) with spin ex- 
change interaction coupled to normal-metal leads can expe- 
rience a quantum phase transition, specifically the singlet- 
triplet transitionr^iiS was recently confirmed by two indepen- 
dent experimentsi^i^ The transition was observed to accom- 
pany a drastic change in the transport mechanism, and it was 
also found that the spin exchange coupling between electrons 
could suppress the Kondo correlation completely or alter its 
physical nature by changing the screening mechanism. The 
influence of such a magnetic interaction on the Josephson cur- 
rent was also studied for a MJJ where the molecule is mod- 
eled by a single-level QD having spin exchange coupling"*^ be- 
tween spins of QD electron and a metal ion. It was predicted 
that the state of the supercurrent can be switched between 
and TT junctions by tuning the magnetic interaction. On the 
other hand, theoretical calculations'* " and experiments^^ have 
shown that the Josephson junction made of a multi-level quan- 
tum dot in the weak-coupling limit can behave as a tt junction 
even when the dot is nonmagnetic without a localized spin and 
vice versa. The studies found out the significant roles of (1) 
the off-diagonal Cooper pair tunneling process^ in which two 
electrons in the pair are transferred via different orbitals in the 
QD and (2) the parity of the QD orbital wave functions" that 
determine the relative sign of the dot-lead couplings. 

In this paper we study the electronic transport through a 
Josephson junction having in it a TLQD with the spin ex- 
change interaction between electrons in two orbital levels. 
Here we focus on the regime where the doubly-occupied QD 
experiences the singlet-triplet transition due to the spin ex- 
change coupling that is tunable by the gate voltage. The 
physical properties of the ground state and the supercurrent- 
phase relation (SPR) through the junction are examined as 
the strengths of the superconductivity and the spin exchange 
coupling are varied. In order to study both of the weak- and 
strong-coupling limits we exploit the numerical renormaliza- 
tion group (NRG) method which can take into account the 
Coulomb interaction in a nonperturbative way. In additions, 
the physical understanding of the numerical outcome is sup- 
plemented by the analytical analysis such as fourth-order per- 
turbation theory and scaling theory. 

Our main findings are summarized as follows: (1) The sign 
of the supercurrent is determined by the competition between 
diagonal and off-diagonal tunneling processes whose strength 
and sign can be controlled by the parity of the orbital wave 
functions and the spin correlation present in the dot. (2) The 
origin and physical property of the TLQD-JJ can be explained 
in terms of the competition between the superconductivity and 
the Kondo correlation found from the normal-lead counterpart 
of the system. For example, the two-stage Kondo effect leads 
to O-TT-O or TT-O-TT doublc transitions with the exchange cou- 
pling or the superconducting gap. (3) When the superconduct- 
ing phase difference between two leads is maximal, the exist- 
ing Kondo correlation is greatly affected. Interestingly, we 
observed that a Kondo-assisted tt -junction can arise if some 
conditions are met. 

This paper is organized as follows: In Sec. II we de- 
scribe the model Hamiltonian of the TLQD-JJ and specify the 
regimes that we are interested in. The weak-coupling limit 



is studied by using the fourth-order perturbation analysis in 
Sec. III. Section IV presents the results of the NRG calcula- 
tions applied to the weak- and strong-coupling limits and the 
phase diagrams of the system with respect to the properties of 
the SPR. In Sec. V we summarize our study. 



II. MODEL 

The TLQD connected to two single-channel s-wave super- 
conducting leads as shown in Fig. [T] is modeled by the two- 
impurity Anderson model: H = Huu + ^ld + Ht, where 

■Hdd = ^(ejnj + Uni^riii) + Unin2 + JSi • S2 (1) 

i 

■Ht = ^ \tii 4^QkM + {h.c 



-kj. 



{h.c. 



(2) 
(3) 



Here c^k^ (di/j,) destroys an electron with energy ek with 
respect to the fermi level and spin /j, on lead £ — L,R (in 
orbital i = 1, 2 on the dot); mw = clk^'^^k^ and n., = 
dl^difj. are occupation operators for the leads and the dot 
orbitals. The Coulomb energy of the strength U is assumed 
to depend on the total number of electrons in the dot. The 
Hund's rule in the dot results in the ferromagnetic exchange 
coupling denoted as J (< 0) between the electron spins = 
■i V ,d], a,.,.'di,,i, where cr are Pauli matrices. The left 
and right leads are assumed to have identical dispersion en- 
ergy Ek and superconducting gap A, while a finite phase dif- 
ference (f) = 0i — is applied between them. The energy- 
independent dot-lead tunneling amplitudes tu hybridize the 
electron states between the dot and the leads, which are well 
characterized by tunneling rates Tif — Trp\tii\'^, where pis the 
density of states of the leads at the Fermi energy. 

Since we are interested in the regime of the singlet-triplet 
transition of an isolated dot, we focus on the parameter region 
in which the dot is doubly occupied. In addition, we consider 
the nondegenerate case with a finite splitting = £2 — ei > 
between two orbitals. Figure [T] displays the energy levels of 
two-electron states of the isolated dot as functions of J: three 
singlet states, |2,0, 0;a) with a = 1,2,3 and three triplet 
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FIG. 1: (color online) (LEFT) Sketch of the TLQD connected to 
two s-wave superconducting leads. (RIGHT) Energy levels EQ,s;a 
of two-electron {Q = 2) states in the TLQD as functions of the 
exchange coupling J. The shaded region is of our interest, where the 
singlet-triplet transition occurs in the ground state. 
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states |2, 1, M), where the states are labeled as \Q, S, M) with 
the charge number Q, the spin S, and the z component of the 
spin M. The singlet states. 



|2,0,0;1) =dl/ijO) 



(4a) 
(4b) 



|2,0,0;2) = 4^4j0) 
|2,0,0;3) = - |0) (4c) 

have the energies, Eq-i — 2ei + U, Eo-2 ~ 2e2 + U, Eq.^ = 
ei + €2 + U — 3 J/4:, respectively, and the triplet states. 



|2,1,1)=4X^|0) 
1 



(5a) 



|2,l,0) = -^(4/i^ + 4/i^)|0) (5b) 



|2,1,-1)=4,4, |0) 



(5c) 



are degenerate with the energy Ei = ei + €2 + U + J/4. 
Due to the finite splitting Se > and the existence of the 
inter-orbital Coulomb interaction, the singlet-triplet transition 
is driven by the competition between the states |2, 0, 0; 1) and 
|2, 1, M) [see Fig.[r|. The bai^e singlet-triplet sphtting is then 
defined by 



Ji(°) = Si-So;i='5e+^. 



(6) 



The external gate voltage Vg can tune the singlet-triplet split- 
ting by affecting the level splitting Ser^^ the exchange cou- 
pling strength or both of them. For simplicity, we as- 
sume that the gate-voltage dependency is implemented only 
through J — J{Vg) and that Se or e.^ are independent of Vg. 
Our simplification can still capture the main physics of the 
system as long as the regime close to the singlet-triplet transi- 
tion is concerned. 

The configuration of the dot-lead coupling is another im- 
portant source that can govern the physics of the system. First, 
the number of the effective channels coupled to the dot can be 
controUed4ii If the condition. 
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tlR 



^2R 



(7) 



is satisfied, the dot-lead coupling matrix has a zero eigenvalue, 
and one of the two channels can be completely decoupled 
from the dot under a proper unitary transformation, result- 
ing in a one-channel problem. This reduction of the effective 
channels then affects the Kondo effect greatly, which will be 
discussed later. Secondly, the phase of the coupling coeffi- 
cients has an influence on the interference and consequently 
on the electron transport through the dotjiii24 Even though no 
magnetic field is applied in our system, the (real-valued) cou- 
pling coefficients can acquire an additional phase tt depending 
on the parity of the orbital wave functions on the doti2^ Two 
distinctive cases can then be conceived: iiLiiRi2Li2R > 
when two orbitals have the same parity and iiLiiRi2Li2R < 
when they have the opposite parities. Taking into account the 
essential impacts of the dot-lead coupling and focusing on the 



consequent qualitative features of system states and electron 
transport, we consider two representative cases in this paper: 



case I : tu 
case II : tu 



t, he 



It 

-- -t 



2R 



It 



(8) 



with 7 < 1. The case I deals with the effective one-channel 
problem with Eq. (|7]l satisfied, while the case II reflects the 
two-channel problem with the negative product of coupling 
coefficients. The effect of asymmetric coupling with respect 
to the orbitals is also examined by setting 7 < 1. Another kind 
of asymmetric junction such as tn^ ^ tm that can happen fre- 
quently in realistic experimental setups like break junctions^— 
is not considered in our study because this asymmetry is ob- 
served to make no qualitative impact on the Josephson current. 

Finally, since we are interested in the low temperature be- 
havior, we concentrate for the most part on the Kondo regime. 
The hybridizations Tu are chosen to be far smaller than the 
particle or hole excitations with respect to the two-electron 
states in order to suppress the resonant tunneling. Specifically, 
throughout our study, we choose ei — —1.61?, 62 — —IAD, 
[/ = Z?, and r = 7rp|ip = 0.051), where the half band width 
D is taken as the unit of energy. Here we have also used the 
particle-hole symmeti-y condition ei + €2 + 3U = 0. 



III. WEAK COUPLING LIMIT: A > Tk 

A. Fourth-Order Perturbation Theory 

First, we consider the weak coupling limit where the super- 
conducting gap A is much larger than the Kondo temperature 
Tk, which will be defined in Sec. |IV] In this case the su- 
percurrent can be calculated via fourth-order perturbation the- 
ory in Ht-^'"' " We apply degenerate perturbation theory that 
takes into account the singlet state |2,0, 0; 1) and the triplet 
states 1 2, 1, A/) simultaneously since they are almost degen- 
erate close to the singlet-triplet transition point of isolated 
dot. Unlike the single-level quantum dot studies^^ where it 
is enough to collect only terms that depend on the phase dif- 
ference (f), on the other hand, one must keep track of all the 
(/)-independent terms in the TLQD study because they con- 
tribute to the renormalization of the singlet-triplet splitting,— 
and the transition point is shifted from its unnormalized posi- 
tion, Ji^"-* = 0. Due to the singlet nature of the Cooper pair, 
there exists no coupling between the singlet and the triplet 
states to any order of the perturbation, and the energy of each 

state is separately shifted: Ea = E^^^ + 5Ea{4>) for a = S,T 



with E^^g'^ = Ea-i and E^^' = Ei. The energy shifts are given 
by 



,(0) 



5E„ 



"^2 r 



3al itl + Pa2 X] lulll 



fiaZ ^Ullt + {Pai - P'ai COS cj)) ^ 



+ {Pab - I3'a5 COS ^)7il71R72L72r] , 



(9) 
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FIG. 2: (color online) Listing of typical virtual tunneling processes 
contributing to SEa- The arrows indicate the direction of the charge 
transfers for processes contributing to (a) fiao, (b) /3ai, (c) Pa2, (d) 



/3a3,(e)/3i'i,and(f) /^W. 



where we have defined = ta /t. Figure |2] shows typi- 
cal virtual hopping processes that contribute to each term in 
Eq. (|9]l. The detailed expressions for the coefficients /3ai can 
be found in the Appendix. 

The effective singlet-triplet splitting then becomes 

Ji(0) = Eri^b) - Es{4>) = Ji^"^ + SJi{(P) (10) 



with 



SJi EE 5Et ~ 5Es = 5Ji^^^ + 5Ji 



,(4) 



(11) 



where JJi'^-' and JJi^"*-' consist of the terms that are propor- 
tional to r and r^, respectively. We find that 6Ji is mostly 
positive in the parameter regime of our interest, favoring the 
singlet formation. The singlet-triplet transition point Jc when 
Ji( J = Jc) = 0, which now becomes 0-dependent, is then 
shifted from its bare value jjf'^ = —ASe = —O.&D to a more 
negative value. It should be noted that the second-order con- 
tribution to 6Ji 



(12) 



is finite in contrast to the previous study of parallel double-dot 
system'*' where the leading contribution is found to be of the 
order of F^. The main difference comes from the character- 
istics of the singlet states in two systems. In the double-dot 
system studied by Choi et al, the two quantum dots, each of 
which is singly occupied, are identical and have no Coulomb 
interaction between them so the lowest-lying singlet state is 
|2, 0, 0; 3), while it is |2, 0, 0; 1) in our system due to the ex- 
istence of the finite splitting and the inter-orbital Coulomb in- 
teraction. The singlet state |2,0, 0;3) has the same charge 
distribution as the triplet states |2, 1, M), so the second-order 
perturbation does not give rise to any additional splitting be- 
tween two states [see Fig. |2](a)]. On the other hand, hav- 
ing |2,0,0; 1) as the lowest-lying singlet states, our system 
can exhibit a rather huge renormalization of the singlet-triplet 



splitting that is of the order of F. This second-order term 
(5 Ji'^-* is numerically found to increase as A is decreased. This 
tendency is opposite to the expectation that the renormaliza- 
tion, which is due to the tunneling of Cooper pairs whose am- 
plitude increases with A, should be weakened as A decreases: 
in other words, limA^o ^Ji — 0- This discrepancy should be 
resolved by the higher-order terms oftheorderof (F/A)" that 
are more involved as A decreases: In fact, the fourth-order 
term SJi-'^^ is observed to become negative for smaller A so 
that the renormalization is diminished. Owing to this opposite 
A-dependencies of (5 Ji'^^ and the other higher-order terms, Ji 
varies non-monotonically with A, which in turns implies that 
the transition point Jc also displays a non-monotonic depen- 
dency on A: see Figs.[3]and|5] 

The supercurrent can be calculated via the derivative of the 
energy with respect to the phase difference 0: 



, '2edEa f . 



(13) 



with 



short 




X ^ 2 2 



3l57lL7lR72L72Rj 
(14) 



As a matter of fact, only the virtual processes in Figs. |2] (e) 
and (f) contribute to the Cooper pair tunneling. The /3^4-term 
[(e)] arises from the diagonal processes where both electrons 
in a Cooper pair travel through either the orbital 1 or 2, while 
the /3^5-term [(f)] from the off-diagonal processes with one 
electron traveling through the orbital 1 and the other traveling 
through the orbital 2. Depending on the order of the sequence 
of electron tunneling and the spin coiTelation of dot electrons, 
the coefficients f3'^^ can acquire a relative minus sign owing to 
Fermi statistics. For the singlet state, one can find that 



/3^4>0 and (3's5<0 



(15) 



The negative sign for /3^g is attributed to the processes with 
one electron traveling through a filled level (orbital 1) and the 
other electron through an empty level (orbital 2). It should be 
noted that it is necessary to take into account the dot electron 
correlation exactly in order to determine the supercurrent sign 
coiTectly. Not all the processes contributing to acquire 
the TT phase: For example, the processes with the intermediate 
state 1 2, 0, 0; 3) acquire no phase at all [see Eq. ( IA29H , while 
their amplitudes are always smaller than those of the other 
processes, and finally /^^^ is negative. For the triplet state. 



P'rpA < and P'rpr^ < 



(16) 



because of the presence of local magnetic moments in both 
orbitals.* Apart from the sign, we have found numerically that 
the off-diagonal processes usually have larger amplitude than 
the diagonal ones: 



for a^S,T. 



(17) 



5 



Hence, when the product 7il7ir72L72R is comparable to 
Si 7ix7iTi, magnitude, the sign of the supercurrent dic- 
tates the sign of the off-diagonal term, or that of the prod- 
uct — 7il7ir72L72R regardless of the spin state: the current 
exhibits the 0(7r)-junction for the negative (positive) product. 
Otherwise, that is, if |7il7ir72L72r| < J2i liLliR^ the diag- 
onal term prevails in determining the sign of the supercurrent 
so that the singlet (triplet) state features the 0(7r)-junction be- 
havior regardless of the sign of the product. 

In the following sections, we identify the system state ac- 
cording to its ground-state spin and the sign of the supercur- 
rent in the A- J plane. We use the labels S and T to denote 
the spin singlet and triplet state, respectively. Since the phase 
transition depends on the superconducting phase difference 
as well, the phase boundaries are located at three different val- 
ues of (j): (red line), tt/2 (green line), and tt (blue line). 
Between = and tt boundaries the system is in the inter- 
mediate state having a stable ground state and a meta-stable 
state. The intermediate states are tagged with a subscript that 
represents the meta-stable state spin. For example, the ground 
state in the state Ts is mostly of the spin triplet, while it is of 
the spin singlet at and near cf) = 0, and the system experiences 
a phase transition from spin doublet to singlet as cf) is varied 
from to TT. The state identification is then supplemented by 
the SPR calculated from Eq. ( fT4l i. classifying whether it is of 
either 0^'-' or tt^''' junctions. Two states with same ground- 
state spin can be distinguished if their SPRs are different and 
the boundary between them will be colored in yellow line. 



B. Case I: 7il7ir72l72R > 
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Figure [3] shows the phase diagrams in the A- J plane in the 
case I for various values of 7. The lower bound of A is set to 
r because the perturbation theory works only when F ^ A. 
For 7=1, the Josephson coupling 



short 



2(r/A)2(2/3;4 + /3;, 



(18) 



is always negative because f3'^^ < and 2|/?^4| < \/3'^^\, and 
the current exhibits the tt -junction behavior, no matter what 
values J and A have [see Fig. [3] (a)]. For 7 <C 1, on the 
other hand, the contribution from the off-diagonal term be- 
comes negligible since 



71L71R72L72R < TiXT^ii = 1 + 7'^, 



(19) 



and the sign of la is governed solely by the /3^4 term. The spin 
singlet state is then of the junction since f3g^ > 0, and the 
singlet-triplet transition accompanies the O-tt transition with 
the intermediate states as shown in Fig. |3] (c). Figure |3] (b) 
shows that for intermediate values of 7, both of the and tt 
junctions can appear in the spin singlet state: the and tt junc- 
tions take place in the regions with larger and smaller values of 
A, respectively. The phase boundary separating two regions 
moves toward the smaller A as 7 is decreased. 



FIG. 3: (color online) Phase diagrams in the A-J plane for the case 
I with 7=1 [(a)], 0.6 [(b)], and 0.1 [(c)]. The phase boundaries are 
located when the ground-state spin is changed at <j) — (red line), 
7r/2 (green line), and tt (blue line). The yellow line separates two 
states with same ground-state spin but different SPRs. Each phase is 
shaded in gray scale according to its SPR: lighter gray for junction 
and darker gray for tt junction. Refer the detailed classification of 
the states to the text. The insets show the SPR for (j) £ [—1", tt] at the 
points indicated by the arrows. Here the value of A is swept from T 
to D, and for A < F the above diagrams are not valid. 



Different strength of the Josephson coupling in the spin sin- 
glet and triplet states gives rise to a discontinuous change in 
the SPR in the intermediate states [see the insets in Fig. O 
and a rapid change in the critical current Ic = \Ia\ across 
the singlet-triplet transition as shown in Fig. |4] The numer- 
ical calculation of the supercurrent finds that the supercur- 
rent is stronger in the spin triplet state than in the spin sin- 
glet state: \It\ > \Is\- In the spin singlet state the diagonal 
and the off-diagonal processes make the opposite contribu- 
tions (/3^4 > > P'g^), resulting in a partial cancellation. 
This is not the case in the spin triplet state in which both pro- 
cesses contribute to the tt junction (/3y4, /3^g < 0). For small 
7, on the other hand, such a cancellation does not make a sig- 
nificant role since the /3^5 term becomes much smaller than 
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FIG. 4: (color online) Critical currents as functions of J in units 
of 1^°^'' for the case I with 7 = 1 [(a)], 0.6 [(b)], and 0.1 [(c)] 
for various values of A: O.ID (top), 0.2D, OAD, 0.6D, and D 
(bottom). The insets display the critical currents as functions of A 
in units of e/h for J — —0.7D (in the spin triplet state) and — D (in 
the spin singlet state). 




FIG. 5: (color online) Phase diagrams in the A-J plane for the case 
II with 7=1 [(a)], 0.72 [(b)], and 0.1 [(c)]. Refer to Fig. [3] for the 
details. 



C. Case II: 71L71R72L72R < 



the /3^4 term, so the critical currents in both spin states be- 
come comparable as can be seen in Fig.|4](c). 

The critical current exhibits a non-monotonic dependence 
on A: see the insets in Fig.|4] In two extreme limits, A ^ F 
and A ^ F, the supercurrent should vanish. The supercur- 
rent, induced by the proximity effect that is proportional to 
A, should vanish in the limit A — 0. In the opposite limit, 
the high energy cost y^A^ + eJ of the quasiparticles created 
during the virtual processes suppresses the current. Conse- 
quently, the critical current has a maximum as a function of 
A. For the intermediate values of 7 when the O-tt transition 
occur in the spin singlet state, the critical current can become 
zero at the transition [see the inset in Fig.|4](b)]. 



The phase diagram and the critical current in the case II 
are shown in Figs. |5] and |6] respectively. In this case one 
of the dot-lead tunneling amplitude changes its sign, making 
the product 7il7ir72L72R negative and accordingly revers- 
ing the sign of the off-diagonal contributions. It then switches 
the junction characteristics from tt to junction for the case 
7=1 when the off-diagonal term prevail over the diagonal 
term: compare Fig. [3] (a) and Fig.|5](a). However, for 7 <C 1 
when the off-diagonal contributions are negligible, the nega- 
tive product does not affect the supercurrent and the phase di- 
agram so much: Fig.[3](c) and Fig.|5](c) are almost identical. 
As a result, for the intermediate values of 7, the additional O-tt 
transition now takes place in the spin triplet state in contrast 
to the case I: compare Fig.[3](b) and Fig.|5](b). Another dif- 
ference from the case I is that the critical current is now much 
larger in the spin singlet state than in the spin triplet state as 
long as 7 is not so small: see Fig.|6] The same argument used 
in the case I applies as well: With the negative product, the 
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FIG. 6: (color online) Critical currents as functions of J in units 
of Jf^°"'* for the case II with 7 = 1 [(a)], 0.72 [(b)], and O.I [(c)] 
for various values of A: O.ID (top), 0.2D, OAD, 0.6D, and D 
(bottom). Refer the guide for the insets to Fig. [5] 



term which is negative cuts down the positive contribu- 
tion from the /3yg term, while in the spin singlet state both of 
two terms contributes to the junction. 

In addition to the properties of the supercurrent, the shape 
of the phase boundaries are also different from those in the 
case I. Figures|5](a) and (b) show that for the moderate values 
of 7 the phase boundaries are much shifted toward the spin 
triplet side, implying that the spin singlet state is being further 
favored. Furthermore, the transition point Jc displays a mono- 
tonic dependence on A and does not approach its bare value 
in the limit A 0, which is contradictory to our expectation 
from the previous weak-coupling argument. The inclination 
to the spin singlet state is accounted for by looking at the /?a5 
term in SEa that is proportional to 7il7ir,72L72R [see the last 
term in Eq. ©J. Numerical calculations observe f3s5 > and 
li^ssl > |/3t5|, which means that with the negative product 
7il7ir72L72R < the spin singlet state is much lowered 
than the spin triplet state [see the Appendix for expressions of 
/3s5 and /^ts]- This term is also observed to make the fourth- 
order splitting term SJi^^^ positive, which is the cause of the 



monotonic behavior of Jc- One may suspect that this favor- 
ing of the spin singlet state in the limit A — ^ is the artifact 
of the fourth-order perturbation close to its limit of validity, 
r/A ^ 1. However, the non-perturbative NRG study in the 
following section finds that the system should be of the spin 
singlet state in the vanishing A limit and that it is attributed to 
the complete screening of dot spins by the two-channel con- 
duction electrons, which will be discussed in details in the 
next section. Considering that the Kondo effect which is re- 
sponsible for the screening cannot be coiTectly captured by 
the perturbation theory, it is quite interesting that it still re- 
flects correct asymptotic behaviors in the limit A — > 0: the 
approaching of Jc to its bare value in the case I and the pre- 
cursor of the disappearance of the spin triplet state in the case 

n. 



IV. STRONG COUPLING LIMIT: A < 

In this section we extend our study to the strong-coupling 
limit by using the NRG method that is known to be suitable 
for the non-perturbative study of the low-temperature proper- 
ties of the impurity system. Even though the standard NRG 
procedure'*^ can be directly applied to the original Hamilto- 
nian, we have introduced a unitary transformation 





= 5 


CLk^ 






= 2? 


dip 


Cbk/i 


CRkp 




dbfj. 


d2p 



(20) 



which makes all the matrix elements of the Hamiltonian real 
in order to boost up the speed of the numerical computation. 
Here the unitary matrices are chosen to be 



V2 



-10/4' 



1 
X 



(21) 



where x — ^ ^"d i in the cases I and II, respectively. Un- 
der the unitary transformation, each part of the Hamiltonian is 
transformed into 



E 

s — a.b 

E 

sk 



^T=E 
ss'k/i 



{esUs + Uus^nsi) + Unafi], + JSa 

~ (-i)^A (cLtcLki+(/i-c-; 

[tss' dl^Cs'k^L + {h.c.)] , 



Sb (22) 
(23) 
(24) 



respectively, where (—1)" = —1 and (—1)'' = 1. Here the 
transformed dot-lead coupling matrix is given by 



ha tbb 



V2t 

V2t 



sm ■ 



cos ■ 



4 4 

• <^ 
7 sm J 7 cos J 



sm ■ 



7 cos ■ 



cosf 
-7 sin I 



case I 



case II 



(25) 



The Wilson's NRG technique'*'''*^ consists of the logarithmic 
discretization of the conduction bands, the mapping onto a 
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semi-infinite chain, and the iterative diagonahzation of the 
properly truncated Hamihonian. Following the standard NRG 
procedures extended to superconducting leads,— we evaluate 
various physical quantities from the recursion relation 



Hn+i = VAHn + ^(/Jat^/.a'+im + {h-c.)) 

Sfl 

- E(-l)^^(/J^+it/sV+u + e^-^-)) (26) 

S 

for N > with the initial Hamiltonian given by 




{d^f s'O^ + [h.C.)) 



{-lY^{fl^fl,^ + {h.c.)) . (27) 



Here the fermion operators fsNti have been introduced as a re- 
sult of the logarithmic discretization of the conduction bands 
and the accompanying tridiagonalization, A is the logarithmic 
discretization parameter (we choose A = 4), and 

1 _ p^-{N+l) 



V(l - A-(2^+l))(l - A-(2^+3)) ' 




(28) 



^ (29) 



withf = (1 + A ^)/2. The original Hamiltonian is recovered 
by 



— = lim £:A-(^-i)/2h^. 

D N^oo 



(30) 



It has been knowni^^ that the logarithmic discretization un- 
derestimates the coupling between the conduction-band elec- 
trons and the dot electrons. In order to avoid this problem, we 
multiply r by a correction factor A\ given by^ii2 

A lnAA + 1 

^^^—A^- ^^^^ 
Within the NRG procedure, the spin of the ground state, 
the occupation (n,), and the spin correlation (Si • S2) can be 
directly calculated from the expectation values of the corre- 
sponding operators. The supercurrent can be also obtained by 
calculating the expectation value 



1='- 
2 



Nr - N, 



(32) 



where Ne = "^''k- I" terms of the fermion operators fsOfj., 
the current expectation value is expressed as 




^short 

with the current matrix defined by 



^ {iss'dl^fs'o>. + {h.C.)) (33) 





1 


tab 


iaa 


Ha ibb 


~ 2t 


tbb 


— tba 



(34) 



The Andreev levels are located from the subgap many-body 
excitations which are identified as the poles of the dot Green's 
functions. 



A. Normal Leads: A = 

In the presence of Coulomb interaction and spin ex- 
change coupling, strong dot-lead coupling can induce nontriv- 
ial many-body correlations that may compete with and even 
suppress superconductivity. A promising candidate of such 
many-body correlations in the QD system is the Kondo effect. 
In order to identify nontrivial correlations in our system and 
to elaborate the analysis of the strongly-coupled Josephson 
junction, it is quite useful to investigate the normal-lead case 
with A = 0. The NRG procedure described above is then 
applied by setting A = and 4> — 0: The latter condition, 
though not being essential, is imposed in order to simplify the 
dot-lead coupling matrix [Eq. (l25Tl1. The normal-lead version 
of our system has been well studied in the literature, so we 
briefly summarize the known theoretical analyses and present 
relevant numerical results in our parameter regime for com- 
parison with the superconducting case. 



1. Case I: Single Channel 

In the case I only the lead-6, that is, the symmetrized 
conduction-band channel is coupled to the dot, with the other 
channel completely detached: see Eq. ( l25T l. The two-level QD 
system attached to a single conduction channel has been well 
studied in the context of the quantum phase transition in a 
vicinity of singlet-triplet degenerate point. In this case the 
system can be mapped onto an exchange-coupled S = 1/2 
Kondo model through a Schrieffer- Wolff transformation:— 

"Hcff = HlD + Ja^a ■ Sbb + Jb^b ' Sfcb + Jj Sa • Sf, . (35) 

Here the Kondo spins Sq and are fictitious QD spins de- 
fined on the basis of the spin singlet state |2, 0, 0; 1) and the 
spin triplet states |2,1,M). Both of the Kondo spins are 
coupled to the localized spin of the conduction channel Sf,;, 
associated to the symmetrized combination of left and right 
leads [see Eq. (l2Tli1: Here we define the localized spins of 
conduction-band electron spins as 



2 5Z 5Z '^skp'''MM"^s'kV- 
kk' /^/^' 



(36) 



The effective spin couplings are, up to linear order in F, 
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FIG. 7: (color online) Energy-resolved transmission coefficient T{uj) 
for a two-level QD coupled to normal leads (case I) with ferromag- 
netic (Ji < Jc) [(a)] and antiferromagnetic (Ji > Jc) [(b)] ex- 
change coupling for various values of J/ D (as annotated). 
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Ja/b — 



Jl 




It 



1 



-£2 - C/ - J/4 ei + 2U - J/4 -ei - [/ - J/4 ea + 2C/ - J/4 



1 



-ei - [/ -e2-U - J/ A -ei - [/ - J/4 



(37a) 
(37b) 



Note that one has Ja ^ Jb as long as 7 7^ 0. 

The ground state of the Kondo Hamihonian, Eq. (|35] | can 
be of the spin singlet or doublet depending on the strength of 
the effective exchange coupling Ji and is known to undergo a 
phase transition at the critical coupling Ji — Jc, which is of 
the Kosterlitz-Thouless-type.^' The ferromagnetic side ( Ji < 
Jc) corresponds to an underscreened S* = 1 Kondo model 
where the conduction electrons screen one of the Kondo spins 
and the remaining 5=1/2 spin then couples ferromagnet- 
ically to the conduction band and becomes asymptotically 
free at low energies. The corresponding Kondo tempera- 
ture Tk{J) decreases with increasing |AJ| = | Ji — Jd [see 
Fig.|7](a)]. On the other hand, on the antiferromagnetic side 
( Ji > Jc), a two-stage Kondo effect takes place for small 
^ J 40,53,54 YixsX, the Kondo effect leads to a screening of one 
of the Kondo spins, Sq with the larger coupling (for exam- 
ple we assume Ja > Jb) which therefore defines the larger 
Kondo temperature Tk- For temperatures lower than Tk, the 
second spin S;, is decoupled from the conduction band. At a 
much lower energy scale (denoted as T^), the effective anti- 
ferromagnetic exchange coupling Ji between Sa and S;, then 
induces the second screening due to the local Fermi liquid that 
is formed on the first spin. is then the Kondo temperature 
of the second spin screened by electrons of a bandwidth ^ Tk 
and density of states 1/{ttTk)'^^ 



rpl 



_Tk_ 
AJ 



(38) 



The second Kondo effect leads to a Fano resonance and makes 
a dip in the energy-resolved transmission coefficient^ 



T(c^) 



ii' fi 



(39) 



where we have introduced the retarded QD Green's functions 

Gij'^(t) = -iQ{t) {{difj,{t),dl^^}) and a coupling matrix 7ij/ 

with 7ii — 1, 7i2 — 721 = 7, and 722 — 7^. As shown in 
Fig.|2](b), the dip becomes widened with increasing A J and 
eventually overrides the Kondo peak until Tj^ « Tk at which 
the Kondo effect completely vanishes. 



2. Case II: Two Channels 

Unless the zero-eigenvalue condition, Eq. (|7| is satisfied, 
the dot is always coupled to both of the two conduction-band 
channels. The low-energy physics of the system is then gov- 
erned by the two-channel two-impurity Kondo model with an 
exchange coupling. Similarly to the case I, the effective spin 



model can be derived via the Schrieffer- Wolff transformation: 

HcS = "HlD + Sq • {Ja Saa + Jb Sfcb + Jabi^ab + Sfca)) 
+ Sft • {Ja Saa + Jb Sfcb — Jabi^ab + Sfca)) 
+ Jl Sa • Sb + 2iJab{Sa X Sfc) • (s„b - Sfc^)- (40) 

The exchange coupling coefficients are given by 



Ja 
Jb 



1 



1 



2r 

\ei + 2U - J/A ' -ei-U - J/ A 
( 1 1 



(41a) 
(41b) 



V-e2-t/-J/4 e2 + 2C/-J/4 
_ \/27r /I 1 \ 

" (v-e2-t/- J/4 + ei+2[/- J/4 ) '^^^''^ 

while one obtains the same expression for J\ as Eq. ( I37bl ). 
Here each of two Kondo spins is coupled to composite local- 
ized spins of conduction-band channels. The effective Hamil- 
tonian, Eq. ( l40t is not convenient for further analysis since 
it contains cross terms (s^fc and Sf,a) that do not conserve the 
channel degrees of freedom. We introduce a unitary transfor- 
mations that diagonalizes the conduction-band spin operator 
in the channel basis that is coupled to Sg for q ~ a.b: 



cos sin 
— sin 1? cos -d 



Cakfi 



(42) 



with i3 = ±itan ^[2 Jaf,/(Ja— Jb)] for q ^ a (upper 
sign) and b (lower sign), respectively. In terms of rotated 

conduction-band spins s^s' = iLkk' L^m' cL^'^MM'Cs'kV', 
the spin exchange terms in the effective Hamiltonian read 



S, • ( Jl Saa + J2 Sbb) (43) 
+ • (J3 Saa + J4 Sbb + J5 (Sab + Sba)) + Jl Sa • Sb, 

where q = a{b) for q — b{a) denotes the index of the Kondo 
spin for which the coupled conduction-band spin operator is 
not diagonalized, and the coefficients are given by 

Jl = J + ^ J sec 2z?, J3 = J + SJ cos 4?9 sec 2^, (44a) 
J2 = J-(5Jsec2z9, J4 = J - (5Jcos4i3sec279, (44b) 
J5 = -2(5Jsin2i5 (44c) 

with J = (Ja + Jb)/2 and SJ = {Ja - Jb)/2. The index 
q is chosen between a and b such that either Ji or J2 is the 
largest among the coefficients. Now the scaling analysis is 
ready with Eq. ( |43] ). Suppose that Ji is the largest one. Upon 
decreasing temperature, the Kondo spin Sq is first screened by 
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FIG. 8: (color online) Scaled NRG eigenenergy flows with the itera- FIG. 9: (color online) Phase shift and linear conductance as functions 
tion N for a two-level QD coupled to normal leads (case II) with (a) of J for a two-level QD coupled to normal leads (case II) with 7 = 1 
7 = 0.5 and J/D = -0.9 and (b) 7 = 1 and J/D = -1.2. [(a)] and 0.5 [(b)]. 



the conduction-band spin Saa, which defines a Kondo temper- 
ature Tk.i- Below this Kondo temperature, the spins Sq and 
Saa are energetically frozen so that the remaining degrees of 
freedom is approximately governed by the exchange coupling, 

J4Sq ■ Sbb- (45) 

The antiferromagnetic coupling will eventually screen out the 
remaining Kondo spin Sq at a lower Kondo temperature Tk.2 
since J4 < Ji. Hence the system undergoes two-stage Kondo 
effects^'' as the temperature goes down: Two Kondo spins 
are screened out one by one since their couplings to relevant 
conduction-band degrees of freedom are different in magni- 
tude. The ground state is of the spin singlet at low tempera- 
tures (T < Tr- 2) due to complete screening, while the partial 
screening in the intermediate temperature Tk.2 < T < Tk,i 
leaves the system in the spin doublet. We have confirmed this 
scaling analysis numerically by examining the RG flow of the 
scaled low-lying eigenenergies in the NRG procedure. Fig- 
ure[8](a) clearly shows that the flow is in the high-temperature 
regime attracted by an unstable fixed point and then goes to 
the stable fixed point in the lower temperature. 

The spin exchange coupling Ji Sa ■ Sb can interrupt the 
Kondo correlation. In fact, we have observed that the two- 
stage Kondo effect ceases to happen if the exchange coupling 
Ji is so antiferromagnetic that Ji > fc^Tx,!. In this regime, 
the Kondo spins are frozen to form a spin singlet by them- 
selves before the conduction-band electrons screen them out. 
Hence, at zero temperature the system undergoes a transition 
between a Kondo state and an antiferromagnetic state as the 
exchange coupling Ji is varied. In contrast to the case 1, how- 
ever, the transition does not involve any change in the spin 
state: The ground state in both states is of the spin singlet. 
Note that the two-stage Kondo effect arises in the feiTomag- 
netic side in the two-channel case while the one in the single- 
channel case happens in the antiferromagnetic side. 

It may be interesting to consider a special case when the two 
Kondo temperatures are equal to each other: Tk,i — Tk,2- 
This can happen when Ja = Jb so that Ji — J4 — Ja + Jab, 
J2 = J3 — Ja — Jab, and J5 — 0, giving rise to the exchange 
Hamiltonian: 

Jl (Sq ■ Saa + Sg • Sfefc^ + J2 (Sq ■ Sbb + Sq • S^a) (46) 
+ Jl Sa • Sb- 



Since Ji > J2, the Kondo spins and Sq are simultane- 
ously screened by the localized spins Saa and s^b, respectively, 
defining a same Kondo temperature Tk. The RG flow in the 
NRG procedure confirms that there exist no unstable fixed 
point and that only one Kondo temperature governs the flow: 
see Fig.[8](b). In our system, the condition Ja = Jb is satisfied 
with 7 = 1 under the particle-hole symmetry condition. 

The transport in the vicinity of the singlet-triplet transition 
of the isolated dot and on both of the antifeiTomagnetic and 
ferromagnetic sides can be analyzed by measuring the linear 
conductance from the NRG calculations. According to the 
Landauer-Biittiker formula in terms of the scattering matrixr^ 
the zero-temperature linear conductance can be expressed in 
terms of the phase shift 6s for each channel: 

G = — sin^ <5, (47) 
n 

with phase difference S = Sa — Sb- We have extracted the 
phase shifts from the energy spectrum in the NRG proce- 
dure by using the fact that the fixed point is described by 
a non-interacting Fermi liquid.''^ On the ferromagnetic side, 
the Kondo screening forms a resonance level on each channel, 
which corresponds to a phase shift 7r/2 in both channels and 
(5 = 0. On the antiferromagnetic side, on the other hand, both 
QD electrons occupy the orbital 1 so that Sa = Si = tt and 
Sb — 0, resulting in S = tt- It implies that the conductance, 
Eq. (|47] i must approach zero on both sides of the singlet-triplet 
transition of the isolated dot while it has a maximum near the 
transition when S = tt/2- Figure |9] shows that the phase dif- 
ference increases rapidly from zero to tt near the transition 
point (J « —0.8D) and that the conductance reaches the uni- 
tary limit when S ~ tt/2- The maximal conductance point is 
shifted with respect to the bare singlet-triplet transition point 
J = —0-8D since the dot-lead correlation favors the spin sin- 
glet state energetically [see Eq. (|37b| i]. As can be seen from 
Fig. |9] the zero-temperature linear conductance does not re- 
flect the presence of two different Kondo scales: the qualita- 
tive feature of the conductance is same for 7 = 1 and 7 < 1. 
Inclusion of Zeeman splitting, finite temperatures, or super- 
conductivity can, however, alter the low-temperature transport 
property dramatically if the relevant energy scale is between 
two Kondo temperatures and one of the Kondo correlation 
with lower Kondo temperature is suppressed. In the next sec- 
tion, we study how it happens in Josephson junctions. 
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B. Superconducting Leads: A / 

Now we investigate the TLQD Josephson junction with 
A ^ 0, considering the single- and two-channel cases sep- 
arately as in the study of normal-lead case. The system state 
is identified as in the study of the weak-coupling regime: re- 
fer to Sec. 1111 Al for the definitions of labels of the system state 
and phase boundaries. Here we introduce a new label D to 
denote the spin-doublet ground state which is missing in the 
weak-coupling regime. 

A series of studies on the single-level QD Josephson 
junction ^^'^' i^-i 6,2 1 ,22 jj^ve already revealed that the compe- 
tition between the Kondo effect and the superconductivity can 
lead to a Kondo-driven phase transition between the Kondo- 
dominant and superconductivity-dominant states and that the 
transition can be driven by tuning the relative strength be- 
tween the Kondo temperature Tk and the superconducting 
gap A. Strong conductivity (A ^ Tk) in the leads en- 
forces the conduction electrons to form Copper pairs by them- 
selves and does not interfere the spin correlation between the 
QD electrons. In the opposite limit (A <C Tk), however, 
the conduction-band electrons in the leads screen out the QD 
spins through spin-flip processes. Hence one can expect that 
as A is decreased the system undergoes a phase transition 
from states that prevail in the weak-coupling limit to other 
states governed by the Kondo effect. Below we find out that 
the superconducting gap introduces an infrared energy cutoff 
to the system and acts like a coherent probe for the Kondo ex- 
citation spectrum. Therefore, we expect that distinct Kondo 
effects in the two cases - cases 1 and 11 - should lead to differ- 
ent phase diagrams in the presence of the superconductivity 
even at zero temperature. 



1. Case I: Single Channel 

Weak Coupling Regime — Figure [TO] shows the phase di- 
agrams in the A-J plane in the case 1. As expected, the 
NRG calculations confirm the results of the perturbation the- 
ory in the weak coupling limit (A/Tk ^ 1). The system 
undergoes the singlet-triplet transition through intermediate 
states, St and Ts, as J is tuned. The phase boundaries be- 
tween them are in perfect agreement with ones found from 
the perturbation theory: compare Figs. [3] and [TOl Not only 
the superconductivity-induced renormalization of the singlet- 
triplet splitting is well reproduced, but also its asymptotic be- 
havior (SJi —i' 0) is correctly predicted in the limit A 
where the perturbation theory breaks down. Note that the 
normal-lead contribution to Ji [see Eq. ( |37bH becomes effec- 
tive for A ^ r, leaving SJi finite. The SPRs calculated from 
Eq. ( |33]) also clearly follow those of the perturbative results: 
compare the insets of Figs. [3] and [10] At 7 = 1, the SPR is of 
the TT -junction regardless of the spin of the ground state, while 
that of the spin singlet state becomes of the 0-junction as 7 is 
decreased. 

Strong Coupling Regime — For smaller A, on the other 
hand, the transition to the spin doublet state takes place, which 
is clearly ascribed to the Kondo effect. The conduction-band 
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FIG. 10: (color online) Phase diagrams in the A-J plane for the case 
1 with 7 = 1 [(a)], 0.4 [(b)], and 0.1 [(c)]. The phase boundaries 
are located when the ground-state spin is changed at = (red 
line), 7r/2 (green line), and n (blue line). The yellow line separates 
two states with same ground-state spin but different SPRs. Refer the 
detailed classification of the states to the text. The insets show the 
SPRs for (fi G [— TT, tt] at the points indicated by the arrows. Here the 
solid lines are guide for eyes. 



electrons in the effective single channel screen out one of the 
two QD spins, leaving the other unscreened. We have ob- 
served that the transition takes place at values of A w Tk 
[see red lines in Fig. [TOl, where Tk is the Kondo tempera- 
ture estimated from the width of the transmission coefficient 



12 



in the normal-lead case [refer to Fig.|7]for 7 = 1]. The phase 
transition is highly dependent on the values of J and exhibits 
asymmetric structure with respect to the sign of Jj. On the fer- 
romagnetic side ( Ji < Jc) the system experiences a transition 
from the spin triplet to the spin doublet state with decreasing 
A, while on the antiferromagnetic side ( Ji > Jc) the singlet- 
doublet-singlet double transition is observed. In the strongly 
antiferromagnetic side ( Ji ^ Jc) there exists no transition at 
all. This J-dependence of the transition originates from the 
fact that the spin exchange coupling affects the Kondo effect 
as discussed in Sec. 11 V ATI First, on the ferromagnetic side 
(Ji < Jc), the Kondo temperature decreases with increas- 
ing I AJ|. It explains the shift of the T-D phase boundary in 
Fig. [TOl toward smaller A with increasing | J|.— On the anti- 
ferromagnetic side ( Ji > Jc), a two-stage Kondo effect with 
two Kondo temperature Tk and takes place for small A J. 
As long as A > T^, the second Kondo effect does not ap- 
pear since the superconducting gap blocks any quasi-particle 
excitation within the gap A. Therefore, for < A < T^, 
one Kondo spin remains unscreened, forming the spin doublet 
state. For A < T^, however. Cooper pairs notice the sup- 
pression of the Kondo resonance level, and their tunneling is 
governed by cotunneling under strong Coulomb interaction, 
restoring the weak-coupling supercurrent in the presence of 
the spin singlet correlation. Hence the observed shape of the 
S-D phase boundary and the reentrant behavior are well ex- 
plained by the fact that the first Kondo temperature Tk de- 
creases with increasing A J as in the ferromagnetic side and 
that the second one decreases with decreasing A J and 
vanishes as A J — )■ 0. 

Finite Phase Difference — The transition boundaries de- 
pend on the phase difference cj), which is responsible for the 
occurrence of the intermediate states. One should note that 
in the presence of finite A and (p both of two conduction- 
band channels are coupled to the QD: All the elements of 
the dot-lead coupling matrix, Eq. ( l25T l become finite. In ad- 
dition, the finite A makes it impossible to decouple one chan- 
nel completely via any unitary transformation. However, the 
effect of the second channel is energetically cut off by the fi- 
nite gap A itself in a sense that the Kondo temperature due 
to the coupling to the second channel is always smaller than 
A. The single-channel argument is thus sufficient to account 
for the phase transitions even at finite cj). The dot-lead cou- 
pling matrix, Eq. ( l25T l then indicates that the couplings tab 
and iftb responsible for the Kondo effect are reduced from 
y/2t to \/2t COS J, which accordingly lowers the Kondo tem- 
perature. The decrease of the Kondo temperature is clearly 
demonstrated in the phase diagrams: The phase boundaries at 
(f) = tt/2 [see green lines in Fig. (TO) are located at smaller 
A than those at = and the regime of double transition is 
shrunken due to the increase of the second Kondo tempera- 
ture [see Eq. (l38Tl1. However, not only the diminished dot- 
lead coupling is responsible for the reduction of the Kondo 
temperature. We have found that a scaling analysis with fi- 
nite A produces exotic terms (like proximity terms) which are 
missing in the normal-lead case. Such terms with finite (j) can 
suppress the Kondo correlation further by twisting the phase 
correlation between two leads. We have observed that at the 



maximally twisted condition, that is, (/) = tt, no Kondo state 
appears at all so that the Kondo state exists only in the inter- 
mediate state. Such a vulnerability of the Kondo effect at a 
maximally twisted phase condition to any finite magnetic per- 
turbation was also observed in the magnetic molecular JJ.— 
The Kondo state may survive the maximally twisted condition 
only if there is no additional magnetic interaction ( Ji = 0) as 
in the single-level QD-JJ. In our system, however, no appear- 
ance of the Kondo state at = tt is observed even along the 
Dt-Ds boundary where the effective splitting is supposed to 
vanish (J k, J^). We attribute it to the fact that the super- 
conductivity shifts the energy levels of the QD spin states and 
induces the 0-dependent singlet-triplet splitting, which seems 
to favor energetically the spin singlet or triplet states over the 
Kondo state. 

SPR andAndreev Levels — Once the Kondo correlation pre- 
vails over the superconductivity, a resonant level is formed 
at the Fermi level, and the Cooper pairs tunnel through the 
Kondo resonant state, resulting in a ballistic 0-junction. To- 
gether with the 0-dependent phase transition, the resonant tun- 
neling makes the curve of the SPR break into three distinct 
segments as soon as the Kondo effect becomes effective, as 
seen in the insets of Fig. [10] The central segment resembles 
that of a ballistic short junctions, while the two surrounding 
segments are parts of the tunneling SPR for the spin singlet or 
triplet states. Since the Kondo state does not occur at = ±7r, 
the SPR keeps the three-segment structure and does not be- 
come of the perfect ballistic junction that was observed in the 
single-level QD-JJ^ 

Figure[TT](a) displays typical variations of the Andreev lev- 
els and supercurrent with J at a fixed value of A/Z? = 10^"^. 
On the Dt-Ds boundary with Ji « Jc [see the middle plots], 
the spin-singlet(red) and triplet(yellow) Andreev levels are de- 
generate in the central segment, while the degeneracy is lifted 
in side segments around (/) = ±7r. Any finite effective singlet- 
triplet splitting, Ji ^ Jc, clearly induces a spin splitting in the 
subgap excitations, lifting the degeneracy, and consequently 
shifts the crossing between the ground state and the lowest ex- 
citation toward = 0, shrinking the central segment. Across 
the crossing, the ground-state spin is changed from 1/2 to 
(1) for Ji > Jc ( Ji < Jc). Similarly, the Andreev levels ex- 
hibit discontinuities in the spectra; for Ji > Jc ( Ji < Jc), 
two outmost Andreev levels with spin 1 (0) in the central seg- 
ment cannot remain as one-electron excitations with respect to 
the spin-0( 1 ) ground state at the transition and are replaced by 
new Andreev levels with spin 1/2. In parallel with the abrupt 
change in the Andreev levels, the SPR shows a discontinu- 
ous sign change (note that / cx —dEa/d(p, as the continuum- 
excitation contribution is negligible), culminating in a transi- 
tion from the to the tt state: accordingly, the states Ds/Dt 
and Sd/Td are of the 0' and tt' states, respectively. As Ji 
grows in magnitude, the central segment shrinks and eventu- 
ally vanishes. The SPR then becomes sinusoidal, which is 
that of a tunnel junction. Once the tunneling junction is fully 
established, stronger singlet-triplet splitting does not lead to 
any qualitative change in the SPR. The observed O-tt transi- 
tion is quite asymmetric with respect to the sign of Ji. First, 
the phase transition in the antiferromagnetic region ( Ji > Jc) 
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FIG. 11: (color online) (a) (LEFT) Andreev levels in units of A 
and (RIGHT) supercurrent / in units of J^''°''* as functions of (f) and 
A/D = 10"^ for various values of J as annotated. Red and yellow 
dots/lines denote the spin singlet and triplet Andreev levels, respec- 
tively, (b) Andreev levels (close to the Fermi level) in units of A as 
functions A / D sXcj> = Q for several values of J as labeled. Here we 
used 7 = 1. 



takes place at Ji < Tk, while the state survives much larger 
ferromagnetic coupling ( Ji < Jc). Second, the antiferromag- 
netic spin splitting gives rise to a double O-tt transition that 
restores the spin singlet correlation at small A. Figure [TT] (b) 
clearly shows that for the antiferromagnetic Ji less than Tk 
(at J — —0.77) the Andreev levels make double crossings 
as A is varied and the spin singlet ground state is restored 
at small A. This is not the case in the ferromagnetic region 
(at J = —0.84) where only one crossing appears nor in the 
strongly antiferromagnetic region (at J = —0.75) where there 
exists no crossing at all. Note that the crossing takes place 
only when the spin-doublet ground state is replaced by the 
spin-singlet or triplet ones: No crossing appears in switching 
between the spin singlet and triplet ground states since each 
of them cannot be a single-particle excitation to the other 

Asymmetric Coupling — The phase diagram is also sensi- 
tive to the asymmetry factor 7 as well. The decrease of 7 
gives rise to the shrink of the intermediate states So and Ds 
in the antiferromagnetic side and the shift of the T-To and 



Td-Dt boundaries toward smaller values of A [compare the 
phase diagrams in Fig.fTO). To understand this behavior, one 
should take a look at the typical dependence of the Kondo tem- 
perature Tk on the lead-dot coupling: Tk oc cxp [— vl/Ftot], 
where Ftot — 2(1 + 7^)r is the total hybridization and A is 
a F-independent constant. As 7 is decreased from 1 to 0, the 
total hybridization decreases from 4F to 2F, which leads to 
the exponential decrease of the Kondo temperature. Hence, 
with the smaller Tk, the transitions from the spin double state 
to the spin singlet or triplet states occur at smaller values of Ji 
or A. The exponential dependence of Tk and T^ on 7 makes 
the intermediate states Sd and Ds almost vanish and hard to 
detect even at 7 w 0.4. 

Another effect of the asymmetry is the appearance of a 
second state in the singlet side in small-A region [see 
Fig. [T0](b)]. This state is not like the one predicted from 
the perturbation theory which works only in the large- A limit. 
In the small-A limit, all the high-order processes should be 
taken into account in order to determine the sign of the super- 
current. Even though the analytical analysis taking all the or- 
ders of the processes is difficult to apply, a rough account can 
be proposed: From the lowest-order terms, Eqs. ( IA28b and 
( IA29I ), one can know that the sign of the supercurrent due to 
each tunneling process is determined based on which interme- 
diate two-electron state appears in the middle of the tunneling 
process. For example, processes with an intermediate spin sin- 
glet state |2, 0, 0; 3) make a positive contribution. In higher- 
order processes, different kinds of intermediate states will ap- 
pear one by one, and a negative contribution can arise only 
for processes with an odd number of the intermediate states 
that invert the sign of the current. Since processes with an 
even number of sign-inverting states should outnumber odd- 
number processes, one can claim that higher-order processes 
are likely to make the positive contribution and consequently 
to favor the state. We observe that this state takes place 
for 7 < 1 and that at 7 = 1 the spin singlet state is purely 
the TT state. We guess that the second state also benefits 
from the weakening of the negative /3gg term with decreasing 
7. Figure [TO] (b) shows that the second state expands to- 
ward large-A region as J is increased. It can be explained by 
the argument that with increasing J the positively contributing 
processes with the intermediate state |2,0,0;3) have larger 
amplitude because the energy cost _Eo;3 — E^^-.i = Se — 3 J/4 
diminishes; see Eq. ( IA29I ). At smaller 7, the two states ex- 
pand further and eventually merge to form a single state, as 
shown in Fig.[T0l(c). 

Critical Current — Figure [12] shows the critical current as 
a function of J/ D for given values of A. For all values of 7, 
the spin triplet state has a larger critical current than the spin 
singlet state, which is consistent with the perturbation results: 
compare Figs. |4] and [12] The critical current is considerably 
boosted up in the Kondo-dominant state and even approaches 
the ballistic short -junction value /^hoit go^^ deep into the 
Kondo state. The critical current has its maximum around the 
singlet-triplet transition point and decreases rapidly in the an- 
tiferromagnetic side. A dip in the critical current is observed 
in the antiferromagnetic side at moderate values of 7. This 
dip happens at the boundary between the it state and the sec- 
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FIG. 14: (color online) Phase diagrams in the A- J plane for the case 
II with 7 = 1 [(a)] and 0.5 [(b)]. Refer to Fig.[TO]for the details. 



FIG. 12: (color online) Critical currents as functions of J/D in units 
of Jf for the case I with 7 = 1 [(a)], 0.4 [(b)], and 0.1 [(c)]. The 
numbers next to the lines represent the values of A/D chosen for 
them. Here the solid lines are guide for eyes. 
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FIG. 13: (color online) (a) Occupation {n2) in the second level of 
the QD and (b) spin correlation (Si -82) as functions of A/D for 
the case I at 7 = 1 and (j> = 0. The numbers next to the lines 
represent the values of J/D chosen for them. Due to the particle- 
hole symmetry, we have (ni) = 2 — (712). Here the solid lines are 
guide for eyes. 



QD occupation and the spin correlation between QD spins. 
For the isolated QD, the perfect spin singlet state dictates 
(ni) = 2, (na) = 0, and (Si • 83) = 0, while (rii) = 1, 
("-2) = 1, and (Si • S2) = 1/4 in the spin triplet state. 
This behavior is well reproduced for large values of A [see 
Fig- ESI ■ The strong superconductivity effectively decouples 
the QD from the leads, and the correlations between QD spins 
are left unpolluted. As A is decreased, on the other hand, 
the dot-lead hybridization interferes the QD spin correlation 
and makes (n^) and (Si • S2) deviate from their bare values. 
Furthermore, they exhibit abrupt changes at phase transitions 
to the spin doublet state; (712) and (Si • S2) exhibit the same 
qualitative dependence on A since the occupation in the or- 
bital 2 is directly related to the formation of the spin triplet 
state. Interestingly, the emergence of the Kondo spin corre- 
lation which screens out one of the spins does not completely 
suppress the local spin correlation between QD spins and even 
help its recovery slightly as A is further decreased. 



ond state where the current vanishes completely. Note that 
the Kondo-driven O-tt transition involves intermediate states 
so that the current does not vanish at the phase boundaries. 

Occupation and Spin Correlation — Finally, we'd like to 
mention about the other ground-state properties such as the 



2. Case 11: Two Channels 

Weak Coupling Regime — Now we consider the two- 
channel case in the presence of the superconductivity. Fig- 
ure [14] displays the phase diagrams in this case. In the large- 
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A limit the phase boundaries and the SPR characteristics co- 
incide well with those obtained from the perturbative theory: 
compare Figs.|5]and[T4] At 7 = 1, the SPR is of the 0-junction 
in both the spin singlet and triplet states, and that of the state 
T switches into the tt -junction as 7 is decreased. 

Strong Coupling Regime — The small-A part of the phase 
diagram, on the other hand, features Kondo-oriented struc- 
tures. First, consider the asymmetric case (7 ^ 1) [see 
Fig. [14] (b)] which exhibits typical phase diagram in two- 
channel case. On the antiferromagnetic side no transition is 
observed to take place. In this regime the antiferromagnetic 
coupling is strong enough that the system is frozen in the spin 
singlet state regardless of value of A. On the ferromagnetic 
side, however, the system experiences a double transition with 
decreasing A: a transition from the spin triplet to the spin 
doublet state is followed by a second transition to the spin sin- 
glet state at smaller A. It reflects the two-stage Kondo effect 
discussed in Sec. II V A2l where two Kondo scales Tk,i and 
Tk,2 are operating. For A > TV,i/2, the conduction-band 
electrons form the Cooper pairs by themselves without affect- 
ing the ferromagnetic correlation formed in the QD. Once A 
is lowered below the larger Kondo temperature Tk.i so that 
Tk 2 < a < Tk,i, one of the QD spins is screened out and 
the other spin that is left unscreened defines the spin doublet 
state. For smaller A < Ta' 1/2^ the remaining QD spin is also 
screened out so that the entire system becomes of the spin sin- 
glet. We have spotted such a double transition at any value of 
(f). Figure [14] (b) shows that the phase difference affects the 
Kondo temperatures in such a way that Tk,i {Tk.2) decreases 
(increases) with increasing (p. The interval Tk.i — Tk,2 is the 
smallest at </) = tt. As in the case I, the modulation of the 
Kondo temperatures can be roughly understood from the de- 
pendence of the dot-lead coupling on [see Eq. (|25])1: The 
phase difference redistributes the amplitude of dot-lead cou- 
plings so that the stronger one get weaker and vice versa as 
is varied from to tt. As a result, with decreasing A the sys- 
tem evolves from the spin triplet state T to the singlet state S 
via all the possible intermediate states including the complete 
spin doublet state D. 

Symmetric Case — In the symmetric case (7 ~ 1), how- 
ever, the phase diagram exhibits quite different features [see 
Fig. [14] (a)]. At = only one direct transition from the 
spin triplet to the spin singlet is observed [see a red line]. It 
is surely due to the disappearance of the two-stage Kondo ef- 
fect at 7 = 1: two QD spins are simultaneously screened at a 
common Kondo temperature Tk [refer to Sec. lIV A"2l . Hence 
the system passes from the spin triplet state to the spin sin- 
glet state either by antiferromagnetic coupling or by Kondo 
correlation. At finite 0, however, the redistribution of the dot- 
lead couplings deviate the Kondo couplings for two QD spins 
from their symmetric point so that one becomes larger and 
the other smaller, recovering the two-stage Kondo effect with 
two Kondo temperatures Tk,i > Tk.2 again. In this case 
Tk.i (Tk.2) increases (decreases) with (f) so that the interval 
Tk.i — Tk.2 reaches its maximum at = tt. As a result no 
complete D state arises as the spin triplet state is changed into 
the spin singlet state. 

One more peculiarity in the symmetric case is a cusp in the 



phase boundary at = tt [see Fig. [T4](a)]. Note that such a 
strange structure is also observed in the asymmetric case [see 
Fig. [14] (b)]. Analytical theory to capture its physical origin 
is not available because the perturbative scaling theory with 
finite A is hard to trace down. Instead, we draw a tentative 
argument from its structural resemblance to that caused by 
the single-channel two-stage Kondo effect. Our NRG calcu- 
lation indicates that at = tt the lower Kondo temperature 
Tk.2 in the ferromagnetic side exponentially decreases as J 
approaches the antiferromagnetic region. As soon as the ef- 
fective singlet-triplet splitting J\ becomes antiferromagnetic, 
the spin exchange coupling JjSa • can then cause the sec- 
ond Kondo effect with the Kondo temperature between 
the spin and the local Fermi liquid formed at the spin as 
explained in Sec. llV ATI Accordingly, we have again a two- 
stage Kondo effect that explains the upward convex shape of 
the unusual phase boundaries very well. The condition that 
it can take place is that Tk.2 < Tj^, that is, the second QD 
spin is screened by the continuous degrees of freedom formed 
at the Kondo resonance level at the first QD spin rather than 
by the conduction-band electrons in leads. In the absence of 
superconductivity this does not arise because Tk.2 is always 
larger than Tj^. However, the NRG calculation suggests that 
this condition is satisfied with finite A and more importantly 
with (f) K, n. 

SPR for the Asymmetric Case — The SPRs shown in the in- 
sets of Fig. [14] clearly reflects the Kondo-driven phase tran- 
sition. The asymmetric case [see Fig. [14] (b)] displays the 
similar evolution of the SPR with respect to the alteration of 
the ground-state spin as observed in the single-channel case. 
The formation of the Kondo-assisted resonant level boosts up 
the tunneling of Cooper pairs and develops the ballistic 0- 
junction. Starting from the tunneling tt -junction in the spin 
triplet state, the SPR then has a shape of the three-segment 
structure with decreasing A. The central ballistic part enlarges 
further with lowering A, and the SPR becomes of the com- 
plete 0-junction as the system enters into the state D. Pass- 
ing through the second transition to the spin singlet state, the 
SPR restores the three-segment form and eventually returns to 
the TT -junction. It is worth noting that while the ballistic 0- 
junction is ascribed to the Kondo resonant tunneling, its sup- 
pression at smaller A is also due to the Kondo effect. The de- 
structive interference between two resonant tunnelings leads 
to a subduing of the tunneling current. The interference is not 
complete (while it is the case in the single-channel case) and 
becomes ineffective as J becomes less negative, resulting in 
an increase of the supercurrent. 

The phase diagram in Fig. [14] (b) shows that there exist 
two spin singlet states S with 0- and tt -junctions, respec- 
tively. The perturbative analysis manifests that the spin sin- 
glet state in the case II pertains to the 0-junction behavior. 
Our NRG calculation shows that it is the case even in the 
small A limit as long as the antiferromagnetic coupling pre- 
vails. However, we found that in the Kondo-dominant singlet 
state the SPR exhibits the tunneling tt -junction behavior as 
shown in Fig. [14] (b), following the SPR of the spin triplet 
state in the large-A regime. The transition between two sin- 
glet state does not involve any intermediate state so that the su- 
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percurrent completely vanishes at the boundary (yellow line) 
between them. 

SPRfor the Symmetric Case — The symmetric case, on the 
other hand, displays exotic features in the SPR as well as in 
the phase diagram. The most intriguing observation is that the 
Kondo-assisted tunneling triggers the ballistic tt -junction be- 
havior. As stated before, the Kondo-driven spin doublet state 
does not form at = but comes into being from (f) = tt 
as A is decreased. The tunneling 0-junction SPR in the spin 
triplet state then transforms into the three-segment one that, 
at this time, has the Kondo-driven ballistic feature at its side 
segments, not in the central part. In addition, the SPR in the 
side segments exhibits the vr-junction behavior: As observed 
in Fig. [TT] the change of the ground-state spin from 1 to 1/2 
makes the Andreev levels cross at the Fermi level, which ac- 
cordingly induces sign change of the supercurrent across the 
crossing point. This ballistic tt -junction can develop only if 
the dot-lead couplings are all comparable and their product 
is negative because only this condition ensures that the spin 
triplet state has the 0-junction and the Kondo transition starts 
at = TT with Tks{<P = tt) > TK,i{(t) = 0). The Kondo- 
assisted vr-junction does not reach its full strength, through, 
and further lowering of A shrinks the side segments, finally 
restoring the 0-junction. Hence, in the small-A limit, the sys- 
tem is of the spin singlet in the state. 

Critical Current — Figure [15] shows the critical current as a 
function of J/D for given values of A. In the large-A limit, 
the critical current is larger in the spin singlet state than in 
the spin triplet state as predicted in the perturbation theory. 
On the other hand, the small-A critical current is observed 



to follow the linear conductance obtained in the normal-lead 
case: compare Figs. |9] and [TS] The critical current exhibits a 
peak exactly where the linear conductance reaches its maxi- 
mum. The peak is not due to the Kondo boosting but orig- 
inates from the competition between the Kondo and the an- 
tiferromagnetic correlations at the singlet-triplet transition as 
revealed in Sec. |IV A 21 However, the Kondo boosting en- 
hances the critical current in the spin doublet state. At a given 
A, the system can be driven into the spin doublet state as J 
is decreased, where the unscreened Kondo correlation opens 
a resonance tunneling. Hence, as can be seen in Fig. [15] (b), 
the small-A critical current features a peak and plateau as J 
is varied. A sharp dip is also identified between them. The 
Kondo-assisted plateau is rather weak in the symmetric case 
[see Fig. [15] (a)] and it disappears in the small A limit since 
no spin doublet state exists in this limit [see Fig. [14] (a)]. To- 
gether with the Kondo-assisted tt -junction, this double-peak 
or peak-plateau structure of critical cuiTent with respect to the 
spin exchange coupling J contrast the two-channel case with 
the single-channel one, so it provides a way to distinguish two 
cases in experiments in which the amplitudes of dot-lead cou- 
plings are not known in priori. 



V. DISCUSSION AND CONCLUSION 

We have investigated the physical properties and the elec- 
tronic transport of two-level quantum dot Josephson junctions 
by focusing on two representative dot-lead coupling configu- 
rations, cases I and II [see Eq. ([SJ]. The fourth-order pertur- 
bation theory applied in the weak coupling limit has revealed 
that the parities of dot orbital wavefunctions, that is, the sign 
of the product of dot-lead tunneling amplitudes, can greatly 
affect the sign of the supercuiTent depending on the spin cor- 
relation present in the dot. The key elements that determine 
the cuiTent characteristics are found to be the existence of a 
localized moment in orbitals which reverses the order of elec- 
trons in Cooper pairs and the competition between diagonal 
and offdiagonal tunneling processes. 

In the strong coupling limit the Kondo correlation competes 
with the superconductivity and the spin exchange coupling 
and the system state is determined by their relative strength. 
We have used the NRG method and the scaling theory based 
on the Schrieffer- Wolff transformed Hamiltonian in order to 
examine the Kondo effect in the normal-lead counterpart of 
our system. The effective single-channel case (case I) exhibits 
three different states - underscreened S = 1 Kondo effect, 
two-stage Kondo effect, and spin singlet state - depending on 
the sign and strength of the spin exchange coupling, while in 
the two-channel case (case II), the system displays two-stage 
Kondo effect and spin singlet state. We have found that the 
numerical results from the NRG method applied to supercon- 
ducting case can be understood in terms of comparison be- 
tween the superconducting gap A and relevant Kondo tem- 
peratures: The Kondo coiTelation found in the normal-lead 
case becomes effective once A becomes smaller than the cor- 
responding Kondo temperature. In this way the superconduct- 
ing gap acts like a coherent energy probe for the excitations of 
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the system otherwise. The competition between the supercon- 
ductivity and the many-body correlations present in the sys- 
tem gives rise to phase transitions which accompany abrupt 
changes in physical properties such as ground-state spin and 
SPR: Among the prominent changes that arise once the Kondo 
effect prevails over superconductivity are the boost-up of the 
supercurrent due to resonant tunneling and the appearance of 
strong 0-junction. Knowing the magnitude of the supercon- 
ducting gap that is under control, therefore, it provides a way 
to measure the magnitude of the important many-body corre- 
lations such as the Kondo temperature or vice versa. 

Beside detecting the system excitations, the superconduc- 
tivity directly alters the system state. In the weak coupling 
limit it renormalizes the singlet-triplet splitting and shifts 
the singlet-triplet transition point which now depends on the 
phase difference as well. In the strong coupling limit the finite 
phase difference influences the interference mechanism and 
suppress or enhance the other many-body effects. Especially, 
at the maximally twisted condition {(j) — it), the Kondo corre- 
lation is completely suppressed in the case I and the Kondo- 
assisted tt -junction is induced in the case II. Knowing that the 
TT -junction in most of cases is usually weak due to its per- 
turbative origin, the latter mechanism opens a way to have a 
strong TT-junction. In the light of quantum computational unit 
which is free of any magnetic control, this tt -junction is more 
promising because the state of the junction can be controlled 
by electric manipulation: the spin exchange coupling can be 
tuned by the gate voltage. 

We expect that our theoretical prediction about phase di- 
agrams and SPRs in TLQD-JJs can be explored in exper- 
iment by using state-of-art fabrication techniques. A re- 
cent experiment'*^ measured electronic transport through Cgo- 
molecular junctions and detected the singlet-triplet transition 
and the accompanying Kondo effects which are consistent 
with existing theories. The same experiment group has ex- 
tended their study to superconducting case^^ where the Al bars 
are attached on top of Au leads coupled to Ceo and the super- 
conductivity is induced in Au leads by the proximity effect. 
The Josephson effect predicted in our paper can be then in- 
vestigated by forming a SQUID, one of which arms contains 
the molecular junction. As implemented in Roch et al.*^ the 
spin exchange coupling can be then controlled by an exter- 
nally applied gate voltage. Or, the superconducting gap can 
be tuned by applying a magnetic field as long as it does not 
suppress the Kondo effect. In this way the switching between 
and TT states with respect to the tuning of J and/or A pre- 
dicted in our calculations could be confirmed in experiment. 
Even without using the SQUID geometry, our theory could 
be tested by measuring a critical current through the junction, 
which should exhibit nontrivial dependence on the spin ex- 
change coupling as discussed above. 
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Appendix A: Energy Shift Coefficients in Perturbation Tlieory 

Here we present the detailed expressions of the coefficients 
f3ai for a ~ S,T are defined in the energy shifts, Eq. ^ in the 
fourth-order perturbation theory. On behalf of readability, we 
introduce the following integrals: 



dx/n 
dx/ir 
dx/ir 



A^AO^ 
A[,{0^ 



dxdy I TT^ 



/(a^,y;0/(2;;Cii)/(a;;Cij) 

dxdyj-K^ 

/(a;,y;0/(^)/(y)/(^;eH)/(x;6,) 

dxdyj-K^ 

dxdy/n^ 

/(^,y;e)/(a^)/(y)/(x;a.)/(y;ei,)^ 



(Al) 
(A2) 
(A3) 
(A4) 
(A5) 
(A6) 

(A7) 



where we have defined f{x) = VT+lc^, f{x; ^) = f{x) + 
^, and f{x,y;^) = f{x) + f{y) + ^ and all the integration 
should be done over the region [— _D/A,_D/A]. All the charge 
excitation energies that appear in the expressions are made 
dimensionless: 

eii = (ei-i?f )/A, ei2 = (e2-i?f )/A, (A8) 
e2i--£;i"VA, 62 = (ii;o;i-i?f )/A, (A9) 

63 = {E^-a - E^^'>)/A, 64 = {Eo;3 - 4°^)/A, (AlO) 



e25-(i?i-i?i"^)/A 



(All) 



for a — S,T, respectively. In terms of the integral macros, the 
coefficients f3si for the singlet energy shift are then expressed 
as 



f3so = -Al 
Psi = AiAii - 



[A[ 



I 12 



(A12) 
(A13) 
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Aiifei) + ^11(61) ^11(62) + ^11(62) 
2 2 
Ps2 = 2AiAn + ^'11(62) + ^11 (62) (A14) 

All (64) + ^11 (64) + An (64) + gn(64) 
2 

5(Aii(65) + i?ii(65) - ^'11(65) - SJifes)) 



/3S3 = 2A1A11 



All(64) + Sii(64) 



(A15) 
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3(^11(65) + (65)) 

2 

/3s4 = 2AiAn - An (61) - ^11(61) - An (62) (A16) 
/3s5 = 2Sn(62) (A17) 



^S4 = -Bn(62) + 2^ 



6 



(A18) 



/3^5 = -2A;i(62) + ^'11(^24) + B[,i^24) (A19) 
-3(A'n (65) + 511(65)), 

where the scaled excitation energies are calculated with Ea = 
Es- The coefficients pTi for the triplet energy shift are given 
by 



/3to = — 



A1 + A2 



{Ai + A2){An + A22) 



(A20) 



(A21) 



An (64) + ^22(64) + 2^12(64) ^ ^'12(64) 



+ 



-Bn(64) + -822(64) 



+ 



An (65) + ^22(65) + -612(65) 
2 

Sn(65) + 5^2(65) 



+ A'i2(65) 



JT2 



(Ai+A2)(An+A22) 



(A22) 



- An (61) + Sn(6i) - A22(6i) + -822(61) 
+ 2(Ai2(6i)- -812(61)) 

- An (62) - -811(62) + A'n(62) + -8^1(62) 

- A22(63) + ^22(63) + A^2(63) + ^22 (63) 
+ Ai2(64) — -812(64) 

-8n(64) + -822(64) - A[^{^2i) - A'22(64) 

2 

- 2Ai2(65) + A'n(65) + A^2(65) + -812(65) 



/3t3 = 



JT4 



_ -811(65) + -822(65) 
2 

(Ai +A2)(An +A22) 
2 

An(6i) + A22(6i) 



(A23) 



-812(61) 



Aii(62) + -8n(62) + A22(63) + -822(63) 



(Ai + A2)(An + A22) 
2 

_ An (64) + A22(64) + 2^12(64) 
4 

3An(65) + 3A22(65) + 2-812(65) 



/3t5 = 2Ai2(6i) + -811(61) + -822(61) 
-811(64) + -822(64) 



(A24) 



(A25) 



+ Ai2(64) + 
-3Ai2(65) 



-811(65) + -822(65) 



Al2(64) i?ll(64)+i?^2(64) ..... 

Pt4 = 2 4 ^ ^ 

3a;2(65) b1i(65) + 5^2(65) 



2 4 

/3t5 = -A'ii(62) - -8li(62) - A^2(63) - ^22(63) 

(A27) 

A'll(64)+A^2(64) 



+ ■ 



+s;2(64) 



3(Al,(g25)+A^2(65)) . 

7; -812(65), 



where the scaled excitation energies are calculated with Ea = 
Et- For better estimation of the sign and the magnitude of the 
supercurrent, fuU expressions for (i'g^i^ and /S^^y^ are given 
below: 



P'sb - 



dxdy/'K^ 



2 



/(a;,y;62)/(a;;ai)/(y;ai)/(a;)/(2/) 61 \J f{x;^ii)f{x) 



dx/n 



^T5 



dxdy/n^ 

mfiy) ^ 



1/2 



3/2 



1 



/(a;,2/;64) /(a;,y,;65)y V/(a;;6i) /(2/;6i)y /(a;,2/;62)[/(a;;6i)] 



1 



1 



f{x,y;C24) /(a;,y;65) 
1 3 



1 



1 



/(.^;6i)/(y;6i) /(^;62)/(y;62) 
2 



dxdy/n'^ 



1 



/(a;, 2/; 64) /(a;,y;65)/ /(a;;6i)/(a;;62) 

2 



1 



1 



1 



1 



1 



/(a;,y;62) V/(a^;6i) /(2/;6i)/ .f{x,y:^23) \f{x;^i2) f{y]£,i2) 

2 



1 1 
+ 



+ 



1 



f{x,y;^24) \f{x;^n) f{y;^i2)J f{x,y;^25) \[f{x;^iiW f{x;^n)f{y;^i2) [/(y;62)] 



+ 



+ 



(A28) 
(A29) 

(A30) 



(A31) 
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